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ABSTRACT 

In  this  paper  we  report  a  new  spectroscopic  ellipsometry  technique  that  overcomes 
much  of  the  ambiguity  associated  with  measuring  an  interface  under  a  film.  For  this 
technique  we  match  the  refractive  index  of  the  overlayer  with  an  immersion  liquid  and  then 
perform  spectroscopic  ellipsometry  at  several  angles  of  incidence.  Essentially,  the 
overlayer  is  optically  (not  physically)  removed,  thereby  rendering  the  ellipsometric 
measurement  sensitive  to  the  interfacial  layer  which  is  often  known  to  be  optically  and 
chemically  different  than  either  substrate  or  film.  The  Si-Si02  interface  resulting  from 
thermal  oxidation  of  Si,  and  the  evolution  of  the  interface  with  annealing  is  studied  using 
the  new  technique. 

INTRODUCTION 

It  is  apparent  that  the  interface  region  between  a  semiconductor  surface  and  a  film, 
particularly  dielectric  films  for  MOSFET  devices,  is  of  crucial  importance. 
Consequentially,  the  interface  region  of  electronically  relevant  films  on  semiconductors  has 
been  widely  studied  by  a  variety  of  techniques  (see  for  example  ref  1  and  refs  1-29  in  ref 
1).  These  techniques  fall  into  a  few  broad  classes  that  will  be  briefly  discussed  with  the 
use  of  Fig.l.  First,  there  are  the  optical  techniques  that  can  access  the  optical  response 
of  the  interface  through  an  optically  transparent  overlaycr,  as  seen  in  Figure  la.  While 
many  successfu1  studies  have  been  done  with  the  optimization  of  this  technique,  particularly 
ellipsometry  studies,  the  key  drawback  is  that  the  optical  signals  from  the  top  interface 
must  somehow  be  taken  care  of  (minimized  or  subtracted  etc.),  in  order  to  extract  purely 
interface  information.  Second,  there  are  techniques  that  can  observe  a  surface  with  the 
overlayer  somehow  removed.  Typically,  the  overlayer  can  be  chemically  etched  or 


Figure  1.  Comparison  of  interface  analysis  methods.  DTTC  QUALITY  INSPECTED  4 

physically  sputtered  or  a  combination.  Figure  lb  indicates  that  the  interface  region  can  be 
altered  by  these  aggressive  processes,  and  that  the  information  obtained  about  a  previous 
interface  contains  an  ambiguity.  There  are  intermediate  techniques  such  as  cross-sectional 
transmission  electron  microscopy,  XTEM,  which  neither  removes  the  overlayer  nor  is 
overly  sensitive  to  extraneous  sources  of  signal.  However,  this  technique  requires 
considerable  sample  preparation,  in  particular  the  sample  thinning,  which  can  alter  the 
interface.  Notwithstanding  the  limitations  of  each  technique,  and  taken  in  total,  these 
techniques  have  contributed  significantly  to  our  understanding  of  the  important  interfaces. 

The  present  paper  deals  with  the  development  and  application  of  a  novel  in-situ 
ellipsometry  technique  that  can  access  the  interface  region  for  an  optically  transparent  film 
on  a  surface14.  The  principle  of  operation  of  the  technique  is  to  immerse  the  film 
covered  substrate  sample  in  a  liquid  that  refractive  index  matches  to  the  overlayer  film  as 
illustrated  in  Figure  lc.  In  this  way,  the  optical  response  consists  only  of  reflections  from 
the  film-substrate  interface  region  thereby  removing  the  above  mentioned  ambiguity.  It 
is  assumed  that  the  interface  region  is  optically  distinct  from  either  the  bulk  film  and 
substrate,  an  assumption  borne  out  in  many  studies.  The  enhanced  sensitivity  of  the 
technique  will  be  demonstrated.  By  way  of  application,  the  Si-SiO*  interface  is  chosen  and 
the  nature  and  evolution  of  this  interface  is  examined  as  a  function  of  thermal  annealing. 


THE  TECHNIQUE 

Index  Matching.  The  ideal  immersion  liquid  should  be  non  reactive  and  transparent 
over  a  wide  spectral  range.  Both  the  refractive  index  and  the  dispersion  must  correspond 
to  the  average  refractive  index  of  the  overiayer.  The  refractive  indices  for  the  immersion 
liquids  were  calculated  using  a  three  term  Cauchy  dispersion  formula,  taking  into  account 
temperature  and  using  literature  values  for  the  various  constants3,4.  Using  liquid  mixtures 
it  is  possible  to  adjust  the  refractive  indexes  of  the  ambient  to  that  for  SiCV  Figure  2 
shows  the  spectral  dependencies  of  the  refractive  index  of  pure  carbon  tetrachloride,  CCI4, 
benzene,  and  mixtures3.  Both  CC14  and  C^H*  are  nonpolar  organic  liquids  that  do 
not  interact  with  SiOj.  Also  in  Figure  2  is  shown  the  calculated  spectral  dependencies  of 
the  refractive  index  of  bulk  and  thin  film  SiQ  calculated  using  a  single  term  Sellmeier 
approximation  after  Jellison3.  Since  n  for  thermally  grown  SiOj  films  is  dependent  on  the 
thickness6,7,  the  use  of  immersion  liquids  with  different  refractive  indexes  corresponding 


to  the  average  refractive  indexes  of  the  films  under  investigation  renders  accuracy 
improvements  possible. 


Figure  2.  The  refractive  index  versus  photon  energy  for  pure  carbon  tetrachloride  (CCU, 
and  mixtures  with  benzene  (QHJ,  along  with  bulk  and  thin-film  SiOj. 


Figure  3.  Model  for  the  substrate-interface-overlayer-ambient  system. 

Analysis.  In  order  to  assess  the  sensitivity  of  the  technique,  the  complicated  structure 
shown  in  Figure  3  and  thought  to  represent  the  Si-thermally  grown  SiO^  system  is 
simplified  to  a  two  film  system  with  interface  and  overlayer  thicknesses  and  Lw, 
respectively.  Thus,  at  this  juncture  the  interface  region  in  Figure  3  is  considered  to  be  a 
homogenous  layer  with  an  effective  dielectric  function,  €W.  The  overlayer  or  bulk  Si02 
film  is  represented  with  an  average  refractive  index  n^,.  The  complex  reflection 
coefficient,  p,  for  the  system  is  given  by: 

p  =  tanYexp(iA)  =  p(A,<M^<^A’VL*Av^  (1) 

where  ¥,A  are  the  ellipsometric  measurables  and  the  terms  in  the  parenthesis  on  the  right 
are  parameters,  some  known  a  priori  and  some  to  be  determined.  The  maximum  change 
in  p,  dp/p,  is  the  condition  of  optimum  sensitivity  in  the  measurables,  A,¥.  It  is  desirable 
to  determine  the  optimum  sensitivity  in  terms  of  the  controllable  parameters:  angle  of 
incidence,  4>,  and  wavelength  of  light,  X.  From  the  analytical  solution  for  optimized 
ellipsometric  measurements  of  interfaces  with  thicknesses,  Lw<  <X/4  in  a  thin  film 


structure*,  it  was  shown  that  the  condition  of  5p/p  divergence  or  maximum  sensitivity  is: 

rtf"*1 *  *  rg”“«p<-2ip)  -  0  (2) 

where  the  r’s  are  the  Fresnel  reflection  coefficients  with  subscripts  corresponding  to  the 
interface  between  the  media  with  numbering  starting  from  0  with  the  ambient  and  0  is 
given  as: 

P  -  -  njsin**)  (3) 

It  should  be  observed  that  the  properties  of  the  interface  do  not  influence  the  best 
sensitivity  conditions. 

If  the  refractive  index  of  the  ambient  is  close  to  the  index  of  the  overlayer,  i.e. 
n<, »  n*,,  then  r01 »  0.  The  optimum  sensitivity  condition  from  eqn(2)  becomes  r12  or,  taking 
into  account  n„  »  nw,  ro2  *  0.  For  the  p  wave  this  condition  is  equal  to  a  minimum  of 
the  ellipsometric  angle,  ¥(</>, X)9. 

We  determine  the  optimal  spectral  range  and  angles  of  incidence  from  a  simulation 
of  ¥(</>, X)  and  A dependencies  using  the  model  in  Figure  3  for  the  Si-SiOj  system 
with  the  following  assumptions:  a.  crystalline  Si  substrate  with  known  dielectric  function9; 
b.  interface  microroughness  (see  Figure  3)  with  effective  height  of  0.2nm  and  composition 
of  50%  c-Si  and  50%  suboxide  SiOx,  with  x  =  0.4  and  the  Bruggeman  effective  medium 
approximation,  BEMA,  was  used  to  calculate  a  dielectric  function  of  the  mixture;  c. 
interface  suboxide  transition  zone  of  0.6nm  composed  of  SiO04;  d.  SiOj  overlayer  with  an 
average  refractive  index  n^CUJ;  e.  an  air  or  pure  CC14  ambient  with  refractive  index  n,. 


Figure  4.  Relative  interface  sensitivity  for  A  and  ¥  for  air  and  CC^  ambients. 

Figure  4  shows  a  plot  of  the  interface  sensitivity  for  both  'fr  and  A  in  terms  of  a  relative 
interface  sensitivity  function  «5A(E)  =  Ao^-A^fE)  for  air  and  CCU  ambients.  Where 
Aq(E)  and  A^/E)  were  calculated  for  the  cases  without  and  with  an  interface  layer, 
respectively.  It  is  seen  that  the  sensitivity  of  A  to  the  presence  of  an  interfacial  layer  is 
increased  by  more  than  an  order  of  magnitude  with  the  parameters  at  optimum  sensitivity, 
but  the  ¥  sensitivity  is  low.  However,  changing  4>  to  78°  for  the  same  system  enhances 
the  ¥  sensitivity  to  the  interfacial  layer  and  this  is  also  shown  in  Figure  4. 

Calculation  Scheme.  In  order  to  find  unknown  parameters  of  the  modeled  interface,  we 
used  the  Marquardt  non-linear  best-fit  algorithm  which  minimizes  the  value  of  the  error 
function: 

<?  -  E  Mfl+pEpP)  (7y  (+„V)  -  T7)J)  (4) 


where  P  is  a  vector  of  N  unknown  interface  parameters  and  Ej  is  the  photon  energy,  and 
the  superscripts  cal  and  exp  refer  to  calculated  and  experimentally  derived  values.  A®*1  and 
are  the  values  obtained  using  the  vector  P,  the  Fresnel  formulas10  and  a  matrix 
algorithm  for  the  multilayer  system  complex  reflection  coefficient.  Optical  parameters  for 
the  immersion  liquid  and  SiOj  overlayer  included  in  the  fitting  procedure  were  calculated 
from  dispersion  equations  mentioned  above.  BEMA  was  used  to  calculate  the  effective 
dielectric  function  for  a  mixture  of  constituents  with  known  optical  properties.  Volume 
fractions  of  the  constituents  are  treated  as  unknown  parameters. 

The  calculation  returns  the  vector  P  when  an  initial  guess,  P°,  is  input,  i.e.  the  program 
returns  a  value  of  P  at  the  local  minimum  of  Q,  near  the  initial  value.  In  addition, 
correlations  between  the  fitting  parameters  are  expressed  in  terms  of  the  correlation  matrix 
of  the  derivatives  3^/3^,,  3A/3pa  and  a  95  %  confidence  limit  was  used  to  calculate  errors 
in  the  fitting  parameters. 

EXPERIMENTAL  PROCEDURES 

A  commercially  available  vertical  ellipsometer  bench  was  modified  to  become 
rotating  analyzer  spectroscopic  ellipsometer  with  the  details  previously  described11  and 
calibrated  according  to  a  published  procedure12  and  the  details  described1.  Calibration 
was  performed  on  a  sample  with  A  *90°  (Si  wafer  with  thermally  grown  Si02  film)  both 
in  air  and  in  the  immersion  cell  in  liquid. 

A  fused  silica  immersion  cell,  as  shown  in  Figure  5,  has  been  designed  for  the 
variable  angle  of  incidence  and  spectroscopic  measurements.  The  most  significant  feature 
of  the  cell  is  that  the  two  optically  fiat  and  annealed  fused  silica  plates,  serving  as  the 
entrance  and  exit  windows,  are  connected  rigidly  to  the  polarizer  and  analyzer  arms,  but 
not  to  the  cell.  The  windows  were  adjusted  to  be  exactly  orthogonal  to  the  incident  light 
at  the  straight-through  position  of  the  ellipsometer  (<£=90°),  in  order  to  avoid  any  deviation 
of  the  incident  light  beam  when  it  passes  through  adjacent  media  with  different  indexes. 
Two  fitted  metal  tubes  in  each  arm  permit  some  lateral  movement  of  the  window  position 
without  a  change  of  the  window  tilt.  The  cell  is  rigidly  attached  to  a  stage  and  connected 
to  the  tubes  using  chemically  inert  flexible  tubing.  A  change  of  <j>  in  the  range  of  67°-90° 
is  possible  while  maintaining  the  window  alignment  precision  during  the  immersion 
measurements. 

In-Situ  Ellipsometer  Immersion  Cell 
Polarizer  Aim  Analyzer  Arm 

'***»•  n*UP» _ tr  Fund  Sflks  C*d 


Figure  5.  Sketch  of  the  variable  angle  of  incidence  immersion  cell. 

Single-crystal  (100)  oriented  2  Ocm  p-type  silicon  wafers  were  cleaned  prior  to 
oxidation  using  a  slightly  modified  RCA  cleaning  procedure13.  The  samples  were 
thermally  oxidized  to  about  20nm  at  800°  C  in  a  fused  silica  tube  furnace  in  clean  dry 
oxygen  which  yielded  MOS  quality  Si02  films  on  Si.  For  the  annealing  studies,  after  a 
particular  oxidation,  one  sample  was  removed  without  annealing  as  a  control,  and  the 
others  were  annealed  in  a  clean  nitrogen  atmosphere  for  annealing  temperatures  and  times 
in  the  range  750-1 100°C  and  1-120  min,  respectively. 

In  order  to  use  the  convergence  routines,  an  estimation  of  the  overlayer  thickness 


is  required.  One  layer  (substrate-overlayer)  and  two  layer  (substrate-interface-overlayer) 
models  have  been  used  to  analyze  the  spectroscopic  ellipsometry  scans  of  the  about  25nm 
SK>2  films  on  Si  samples  in  air.  From  these  measurements  in  air,  the  difference  in  the 
minimum  error  function,  Q,  for  one  and  two  layer  models  was  negligible,  which  indicates 
low  interface  sensitivity.  However,  these  measurements  yield  a  good  estimate  of  the 
overlayer  thickness  to  a  precision  of  better  than  0.5%.  The  refractive  index  of  the  film 
was  calculated  from  dispersion  relations.  Immersion  measurements  were  performed  in 
pure  CCI4  at  20°C.  Carbon  tetrachloride  becomes  nontransparent  at  energies  higher  than 
— 4  eV  which  then  defines  the  upper  limit  of  the  spectral  range  to  be  4eV.  The  precision 
of  rotating  analyzer  ellipsometry  without  using  an  achromatic  compensator  falls 
significantly  when  A  approaches  180°  or  0°  w.  To  avoid  this  situation,  we  have  used  2.5 
eV  as  the  lower  energy  limit  where  A  is  more  than  20°.  Therefore,  A  is  restricted  to 
20°<A<160°  in  the  spectral  range  2.5-4. 0  eV  for  the  measurements  in  the  immersion 
liquid  at  the  angles  of  incidence  <t>  —  70°-80°. 

RESULTS:  Evolution  of  the  Si-SiOj  Interface  During  Annealing. 

Model  Independent  Results.  A  series  of  annealing  experiments  were  performed  on  SiOj 
covered  Si  samples,  and  the  results  are  first  displayed  without  recourse  to  a  model,  so  that 
an  objective  view  of  the  interface  evolution  can  be  obtained.  Figures  6  and  7  display 
experimental  results  in  terms  of  the  interface  parameter  SA^XT^tjJ  defined  as: 

-  H’SVj'J  <5> 

where  A^fT.,  O  is  the  experimental  ellipsometric  angle  A  at  a  specific  annealing 
temperature  and  time,  A„  is  the  ellipsometric  angle  for  a  non-annealed  sample  and  the  term 
SA^CT,.,  t„J  is  the  overlayer  relaxation  correction.  This  correction  term  represents  the 
difference  in  A  for  the  single-film  (without  interface)  system  of  a  non-annealed  and 


Figure  6.  The  experimental  dependence  of  the  effective  interface  parameter,  A**,  on 
annealing  time  at  a  number  of  annealing  temperatures.  A  simulated  dependency  is  shown 
with  the  dashed  curve. 

annealed  overlayer  with  the  refractive  index  calculated  from  a  consideration  of  stress 
relaxation6,13,16,17;  and  represents  only  a  few  percent  of  8A&.  The  chosen  photon 
energy  E=3.18  eV  is  in  the  range  of  the  maximum  interface  and  the  minimum  overlayer 
sensitivity.  Figure  6  shows  the  values  for  5A w  measured  at  3.18  eV  versus  anneal  time 
at  a  number  of  temperatures,  and  versus  anneal  temperature  at  two  times  in  Figure  7. 
Other  information  relating  to  modelling  is  also  included  in  Figure  6,  but  discussion  of  these 
items  will  be  deferred.  It  is  seen  in  Figure  6  that  at  all  annealing  temperatures  two 


temporal  regions  of  behavior  are  present.  Initially,  the  anneals  yield  a  relatively  fast 


Figure  7.  The  experimental  dependence  of  the  effective  interface  parameter,  A^,  on 
annealing  temperature  for  two  annealing  times. 

increase  in  SA^  that  slows  after  5  min  and  then  nearly  saturates.  The  data  in  Figure  7  is 
from  Figure  6  at  two  anneal  times,  one  before  and  one  near  saturation,  and  this  plot  versus 
anneal  temperature  reveals  a  distinct  break  in  the  900P-970°C  range,  which  corresponds  to 
the  viscous  flow  range,  i.e.  a  temperature  above  which  viscous  flow  of  the  oxide  is  fast18. 
Thus,  we  consider  that  the  viscous  relaxation  dominates  at  the  higher  anneal  temperatures, 
but  at  lower  temperatures  where  a  large  part  of  the  change  in  occurs,  we  need  to 
consider  other  possible  mechanisms. 

If  the  interface  region  is  treated  as  a  single  homogenous  film,  then  the  extent  of  the 
interface  is  observed  to  decrease  with  both  annealing  time  and  temperature  as  evidenced 
by  the  increasing  A^.  With  both  a  decrease  in  the  interface  region  for  short  and  long 
times  and  low  and  high  temperatures,  the  model  that  is  chosen  for  the  different  modes  of 
behavior  must  show  the  observed  functionality. 

An  Interface  Model.  Previous  models  for  the  Si-SiOj  interface  have  agreed  that  the 
interfacial  region  is  structurally  and  chemically  different  from  either  Si  or  SiOj. 
Components  of  the  interfacial  region  include  roughness  and  a  suboxide,  SiOx,  as  is  shown 
in  Figure  3.  The  parameters  that  quantify  this  working  model,  and  then  the  experimental 
data  A^CE.^.T^.O,  ir“p(E,<£,T,0,t-)  are  reduced  to  yield  values  for  the  model  parameters 
which  are  compared  with  other  studies  about  the  interface  region.  We  describe  the 
roughness  by  crystalline  silicon  protrusions  as  hemispheres  with  an  average  radius  R, 
which  form  a  hexagonal  network  with  an  average  distance  D  between  centers.  The 
protrusions  and  the  region  between  them  are  covered  by  a  layer  of  suboxide,  SiOx,  with 
0<x<2,  and  with  an  average  thickness  L*.  The  effective  thickness  of  the  complex 
transition  layer  is  given  as: 

LW‘R  +  L»,  (6) 

and  an  effective  dielectric  function  e*,,  which  represents  a  mixture  of  crystalline  silicon  c- 
Si,  silicon  suboxide  SiOx  and  the  Si02  overlayer  and  written  as: 

eiV  “  eii^ee-Si*fe-5P  €SIOg'faO*  ^SlOjfsio) 

was  calculated  using  the  BEMA,  where  the  dielectric  properties,  e,  and  relative  volume 
fractions,  f,  of  all  of  the  interfacial  layer  constituents  are  known  a  priori.  The  dielectric 
function  of  SiOx  was  calculated1,2  using  the  BEMA,  and  by  considering  that  SiOx  is  a 
mixture  of  a-Si,  and  SiO*1’. 


In  order  to  model  the  evolution  of  the  interface  during  annealing,  we  use  a  power 
law  to  describe  the  reduction  of  both  tne  protrusions  and  the  chemical  transition  layer  with 
the  powers  p  and  g  respectively  are  taken  to  be  p=g=0.5  which  implies  a  diffusion  model 
and  is  justified  below. 

The  minimization  of  the  error  function  (eqn.(4))  for  the  sets  of  experimental  data  for 
the  non-annealed  wafer  at  and  75°  gives  the  average  distance  between  the  centers 

of  protrusions  D  =  44 ±4  A,  initial  radius  of  the  protrusions  R°  =  9.8±0.3  A  and  initial 
thickness  of  suboxide,  SiO,  L*  =  3.4±0.2  A.  These  results  are  in  agreement  with  the 
interface  geometry  sensitive  TEM  study  of  the  SiO^-Si,  which  shows  that  distances  between 
protrusions  at  non-annealed  interface  are  distributed  in  the  range  40-50  A  and  heights  are 
9-15  A20.  Also,  photoelectron  spectroscopy  revealed  a  chemical  transition  layer  with  a 
thickness  of  2.4-4  A21. 

A  simulation  of  the  kinetic  dependencies,  SAj^T.,  t*J,  with  p=g=0.5  yields  close 
agreement  (dashed  lines  in  Figure  6)  with  the  experimental  plots  of  SA^t^). 


Figure  8.  The  dependence  of  the  effective  interface  thickness  on  annealing  time  at  a 
number  of  annealing  temperatures. 

Figure  8  displays  the  decrease  of  the  extent  of  the  interface,  L^,  i.e.  effective 
thickness  from  eqn.(6)  with  a  minimum  thickness  of  about  3  A  realized  at  the  highest 
anneal  temperatures.  The  interface  effective  refractive  in^ex  vs  t^  at  several  TM  is  shown 
in  Figure  9. 

Initially  both  the  index  and  thickness  of  the  interface  layer  decrease  due  to  the  shrinking 
of  the  silicon  protrusions.  This  is  followed  by  the  slower  decrease  of  the  SiO*  transition 
layer  which  becomes  dominant  aiter  considerable  reduction  of  the  height  of  the  protrusions. 
The  sharp  increase  in  the  index  seen  in  Figure  9  for  the  1 100°C  annealed  sample  is 
indicative  of  the  index  returning  quickly  to  the  SiO  value  as  the  fraction  of  the  Si  as 
protrusions  in  the  interface  layer  rapidly  goes  to  zero.  The  same  effect  would  occur  for 
the  other  anneal  temperatures  but  more  slowly.  A  slower  rise  in  index  is  seen  for  the 
1050°C  annealed  sample  after  10  min  and  the  time  is  too  short  to  see  the  rise  for  the  lower 
temperature  anneals.  The  overall  shape  of  the  results  in  Figure  9  is  typical  of  and  dictated 
by  the  BEMA  model  used  to  interpret  the  data. 

The  results  show  two  distinct  modes  for  the  evolution  of  the  Si-SiQ  interface  upon 
annealing.  The  low  temperature  mode  corresponds  to  a  reduction  in  the  suboxide  and 
protrusions.  The  coefficients,  p  and  g  that  scale  this  evolution  yield  an  excellent  fit  to  the 
data  when  a  value  of  0.5  is  used  in  modelling.  This  suggests  that  Si  atom  diffusion  is 
operative  for  the  smoothening  reaction  and  inward  trace  oxidant  diffusion  for  the  suboxide 
disappearance.  At  high  temperatures  these  reactions  occur  but  the  dominant  mode  is 
viscous  relaxation  of  the  oxide  overlayer.  These  model  dependent  assertions  are  consistent 


with  other  studies  and  are  shown  to  be  accessible  by  the  novel  ellipsometric  technique 
herein  described. 


Figure  9.  The  dependency  of  the  effective  interface  refractive  index  on  annealing  time  at 
a  number  of  annealing  temperatures. 
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